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ABSTRACT 

A difficult and long-standing problem in mathematical physics concerns the determi- 
nation of the value of f{oo) from the asymptotic series for f{x) about x — 0. In the past the 
approach has been to convert the asymptotic series to a sequence of Pade approximants 
{P^{x)} and then to evaluate these approximants at a; = oo. Unfortunately, for most 
physical applications the sequence {P^(oo)} is slowly converging and does not usually give 
very accurate results. In this paper we report the results of extensive numerical studies for 
a large class of functions f{x) associated with strong-coupling lattice approximations. We 
conjecture that for large n, P"(oo) ~/(oo) -f- B/lnn. A numerical fit to this asymptotic 
behavior gives an accurate extrapolation to the value of f{oo). 
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There are many examples in mathematical physics where it is necessary to calculate 
the value of a function f{x) at a: = cxo but where f{x) can only be determined perturbatively 
for small x.^''^'^ The result of such a perturbative calculation yields the first few terms of 
the asymptotic series for f{x): 

oo 

fix) ^ X'^Y. ^r.x'^ ■ (1) 
n=0 

The problem is then to use this series to determine f{oo) assuming, of course, that f{oo) 
exists. 

A natural approach^'^'^'^'^ to solving this problem is to raise both sides of (1) to the 
power 1/q!, 

oo 

f{xf/- = xJ2 <^nX- ' (2) 

and then to convert the right side of (2) to a finite sequence of Pade approximants, P^{x), 
n= 1, 2, 3, . . . , AT, where 2A'" — 1 is the highest order of the perturbative calculation used to 
obtain (1). Since P^{oo) exists, we compute the sequence of Pade extrapolants [P^(oo)]" 
with the hope that as n increases, this sequence rapidly approaches /(oo). 

While this extrapolation procedure works well in some cases one encounters an im- 
portant class of functions for which this procedure is rather ineffective. For this class of 
functions the coefficients bn in (1) have two characteristic properties: 

(i) bn ^ n\ as n ^ oo; 

(a) the bn exhibit a doubly-alternating sign pattern +, +, — , — , +,-!-,—,—,... . 
Asymptotic series whose coefficients exhibit properties (i) and (u) arise in strong- coupling 
lattice expansions in quantum field theory.^'^'^'^'^ The problem of determining the value of 
/(oo) corresponds to extrapolating to the continuum limit of the lattice theory. For such 
functions, experience shows that while [P"(oo)]" often comes close to /(oo), the sequence 
of Pade extrapolants appears to level off at some value that differs from /(oo) by several 
percent. This raises the question of whether the sequence actually converges to /(oo) or 
to some other limit near /(oo). 

Here is a simple example that illustrates this problem: Consider the series 

oo 

f{x) = J2 ^nX^ ' (3) 
n=0 

where 

a2n = (-1)" (2n)! , as^+i = (-1)" 2 (2n)! . (4) 

We can evaluate this series exactly using Borel summation. The result is the integral 
representation 

oo 

m = (1 + 2x) J dt , (5) 
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from which we determine that 

/(OO) = TT . (6) 

Let us try to obtain this result numerically by converting the series in (3) to a sequence 
of Pade extrapolants. On Fig. 1 we plot the value of P^(oo) as a function of l/n for 
n = 1, 2, 3, . . . , 20. Observe that the Pade approximants form a zigzag sequence^ that 
seems to tend to a limiting value on the y-axis (corresponding to n — oo). If we extrapolate 
the sequence to n = oo using a straight line, the sequence appears to have a limit that 
is too small (roughly 10% smaller than n). However, if we include many more terms 
(n = 1, 2, 3, ... , 300) in the sequence (see Fig. 2), we can see that the sequence actually 
lies on a curve that becomes steeper as n ^ oo. Figure 2 indicates that the asymptotic 
behavior of P^{oo) for large n cannot have a straight- line form tt + B/n. 

Because the sequence in Fig. 2 lies on a curve and not on a straight line it is difficult 
to predict its limiting value. A detailed numerical fit of the form 

p:{oo) - 

does not work for any value of C. However, a more general numerical fit of the form 

B 



n*-" (Inn)^ 

gives C = ± 0.1 and i:' = l ± 0.1. Thus, we believe that 

P,"(oo) ~ /(oo) + ^ . (7) 

In fact, a fit of the form P^{oo) — A ^ B / \n.n gives A = 7r with a relative error of less 
than 1%, which is an order-of-magnitude improvement over the result of extrapolating the 
curve in Fig. 1. In Fig. 3 we plot the Pade extrapolants shown on Fig. 2 as a function of 
l/ln(20n).^° Observe that the zigzag curve now lies along a straight line that appears to 
approach the value tt. 

Next, we consider the more complicated case of the Stirling series expansion for the 
factorial function: 



-1 M 1 1 2 139 3 571 4 

fix) ~ a; 2 1-1 X -\ X X x 

''^ ' ^12 288 51,840 2,488,320 

I 163-879 5,246,819 

209,018,880 75,246,796,800 " " ^ ' ' 



(8) 



where 

f{x) = ^e^x^ (i)! . (9) 
V27r X 

Note that the series in (8) exhibits properties (i) and (ii).^^ 
Since 0! = 1 we see that 

/(oo) = = 0.398942... . (10) 

V 27r 
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We attempt to calculate /(oo) by raising (8) to the power —2, computing -P^(x), and 
forming the sequence of extrapolants [P"(oo)]~^/^. As in the previous example if we plot 
these extrapolants as a function of 1/n for n = 1, 2, 3, . . . , 90 (see Fig. 4) we see that the 
extrapolants do not lie along a straight line. From just the first 20 extrapolants we would 
predict the value of f{oo) to be about 0.42, which is 5% high. However, it is clear that 
the extrapolants follow a curve which becomes increasingly steep as n — > oo. If we replot 
the data as a function of l/ln(20r2,) (see Fig. 5) we see that the same extrapolants lie on 
a straight line that appears to intersect the y-axis at precisely the correct result in (10).^° 
A numerical fit of the form in (7) gives /(oo) correct to about 1%. 

Our third example is taken from a quantum field theory calculation. Consider a g(f)^^ 
self-interacting scalar quantum field theory in £)-dimensional Euclidean space. In recent 
papers^^'^^ we showed how to express the free energy for such a theory as a series in 
powers of D in the limit of strong bare coupling, (7 ^ oo. The coefficient of D in this 
dimensional expansion can be expressed as the logarithm of an infinite series in powers of 
X like that in (1); we denote this infinite series f{x).^'^ For all K the coefficients of the series 
f{x) grow like n! and have a doubly-alternating sign pattern. The expansion parameter, 
x = [ga'^^~-^^^~^^~^/^ J is a dimensionless combination of the coupling constant g and 
the lattice spacing a, which is small for fixed a and large g. But in the continuum limit 
for which a — 0, we have a; — ^ 00. Thus, we are confronted with the problem addressed in 
this paper, namely that of computing /(oo). 

As a special case, we consider the free theory corresponding to K = 1. For this case 
there is a closed- form expression for f{x): 

f{x) = X exp|2^ dte-Mn[e-^%(a;t)]| ,. (11) 

where Iq is the modified Bessel function of order 0. In the continuum limit, x — > 00, we 
have 



/(oo) = ^ = 0.28347... . (12) 
The first few terms of the asymptotic series representation for f{x) in powers of x are 

/ ^ o 10 o 29 4 11 . 391 fi 2389 7 5303 o 

f(x) = x[l-2x + 3x^ x^ + —x^ x^ + x^ x^ - ——x^ 

' \ 3 12 10 180 630 448 /-,o^ 

2602051 9 159662191 651255947 ^ _ 435388434359 12 \ 

^ 90720 ^ ^ 907200 ^ 1663200 ^ 119750400 ^ +■■■) 

If we convert the first 250 terms in this series to a sequence of Fade extrapolants, we 
obtain results similar to those of the two examples above. On Fig. 6 we plot this sequence 
of extrapolants as a function of 1/n and observe that the extrapolants lie along a curve 
that is difficult to extrapolate to its value at n = oo. However, on Fig. 7 we plot these same 
extrapolants versus 1 / ln(4n) and obtain a straight line whose limit is clearly very close to 
the exact answer given in (12).-*^^ 

On the basis of our numerical studies of these examples we conjecture that the asymp- 
totic behavior in (7) of the nth Fade extrapolant is generic if the coefficients exhibit prop- 
erties (z) and (ii). We may express the Fade extrapolant P"(oo) constructed from the 
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asymptotic series x X^^o ^n^^ cis a ratio of two determinants, an (n + 1) x (n + 1) deter- 
minant divided by an n x n determinant.^^ Thus, a simple way to state this conjecture is 
in terms of determinants of large-dimensional matrices: 



C(oo) = 
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/(oo) + 



B 

Inn 



(14) 
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FIGURE CAPTIONS 



Figure 1. The first twenty Fade extrapolants P^{oo) for the series whose coefficients are 
given in (4) . The approximants are plotted as a function of 1/n such that the vertical axis 
corresponds to n = oo. Successive Fade extrapolants are connected by line segments to 
form a zigzag line. The vertices of the zigzag line are the values of the approximants. One 
would think that it is easy to extrapolate by eye to the limit of the sequence (the point 
where it intersects the vertical axis). However, the zigzag line of extrapolants appears 
destined to intersect the vertical axis well below the exact answer, tt, which is indicated 
by the horizontal line. 

Figure 2. Same as in Fig. 1 except that we have plotted the first 300 approximants. 
Having plotted many more approximants than in Fig. 1, we believe that the zigzag line of 
approximants is curving upward; it could well intersect the vertical axis at tt. It is difficult 
to guess where such a curve might cross the vertical axis. 

Figure 3. The same approximants as in Fig. 2 now plotted versus l/ln(20n). The zigzag 
line of approximants seems to be tending to tt, the correct value of /(oo). 

Figure 4. The first 90 Fade extrapolants constructed from the Stirling series (8) for the 
factorial function. As in Fig. 2, the zigzag line of extrapolants, when plotted versus 1/n, 
begins to curve as n gets large. For small values of n the zigzag line appears to be headed 
for an intersection with the vertical axis that is too large compared with the correct answer 
l/-\/27r = 0.398942 ... . For larger values of n, the curvature of the zigzag line makes it 
too difficult to predict where the line will cross the vertical axis. 

Figure 5. The same approximants as in Fig. 4 now plotted versus 1/ ln(20n). As in Fig. 3, 
the zigzag line of approximants seems to be tending to the correct value of f{oo), which 
in this case is l/A/27r== 0.398942. Apparently, one can use the Stirling series for (1/a;)!, 
which is valid clS X — 0, to calculate 0! = 1. 

Figure 6. The first 125 Fade extrapolants constructed from the asymptotic series (13) taken 
from a field-theoretic calculation. As in Figs. 2 and 4, the zigzag line of extrapolants, when 
plotted versus 1/n, begins to curve as n gets large. For small values of n the zigzag line 
appears to be headed for an intersection with the vertical axis that is too large compared 
with the correct answer e'^ / (27r) = 0.28347 .... For larger values of n, the curvature of the 
zigzag line makes it too difficult to predict where the line crosses the vertical axis. 

Figure 7. The same approximants as in Fig. 6 now plotted versus l/ln(4n). As in Figs. 3 
and 5, the zigzag line of approximants seems to be tending to the correct value of /(oo), 
which in this case is 3.52775 ... . 
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